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ABSTRACT 

The orbital parameters of extra-solar planets have a significant impact on the probability that the 
planet will transit the host star. This was recently demonstrated by the transit detection of HD 17156b 
whose favourable eccentricity and argument of periastron dramatically increased its transit likelihood. 
We present a study which provides a quantitative analysis of how these two orbital parameters affect 
the geometric transit probability as a function of period. Further, we apply these results to known 
radial velocity planets and show that there are unexpectedly high transit probabilities for planets at 
relatively long periods. For a photometric monitoring campaign which aims to determine if the planet 
indeed transits, we calculate the expected transiting planet yield and the significance of a potential 
null result, as well as the subsequent constraints that may be applied to orbital parameters. 
Subject headings: planetary systems - techniques: photometric 



1. INTRODUCTION 

With the number of known extra-solar planets exceed- 
ing 300, statistical interpretations of the distribution of 
orbital parameters are becoming increasingly significant. 
These parameter distributions help us unlock the mys- 
teries surrounding the planet formation process to which 
many challanges have been presented, not the least of 
which contains the mecha nisms that drive planet a ry mi - 
gration (Armitage 200^. iFord. Ouinn. fc VerasI (|2008l ) 
showed that transit light curves in particular can be used 
to characterize orbital eccentricities and hence give fur- 
ther insight into the global eccentricity distribution. 

In terms of the sheer number of transit light curves, 
the major contributors have been the shallow wide- 
field surveys such as the Tran satlantic Exoplanet Sur- 
yey (TrES) (iMandushev et all [2007i). the XO project 
pohns-KruU et al.l |2008[ ). the Hungarian A utomated 
Telescope Network (HATNetl (IPlIet al.l l2008D. and Su- 
perWASP (j Anderson et all l2008f l. hTaddition, there 
have been at least five cases in which planetary transits 
were detected through photometric follow-up of planets 
already known via their radial veloci ty (RV) discoveries. 
These five planets are H P 209458b (ICharbonneau et a.1.1 
[2000t iHenrv et al.l[2000h. HD 149026b (iSato et al.ll2005n. 
HD 189733b (BouchiLet ah 2005), GJ 436b (Gillon'eraLl 
|2007( ). and HD 17156b (Barbieri ct al. 2007 ). The case of 
HD i7156b is of particular interest since it is a 21.2 day 
period planet which happens to have a large eccentricity 
(e ~ 0.67) and an argument of periastron which places 
the periapsis of its orbit in the direction toward the ob- 
server and close to parallel to the line of sight, resulting 
in an increased transit probability. 

Conversely, the dominant sources of RV planet discov- 
eries have been the California & Carnegie Planet Search 
(|Marcv et al.lll997l ) and the Hig h Accuracy Radi al veloc- 
ity Planet Searcher (HARPS) (,Pepe et al.ll2004[ ) teams. 
However, in t he near future we can expect to see larger- 
scale surve ys (iKane. Schneider. fc~G3[200l and new in- 
struments (|Li et al.ll2008[ ) which wHl increase both the 
number and diversity of known planets. There have 
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been suggestions regarding the strategy for photomet- 
ric follow-up of these rad ial velocity planets at predicted 
transit times (|Kand[2007l) and the instruments that could 
be used for such surveys ()L6pez-Moralesll2006l ) . Some at- 
tempts have been made to detec t these possible t ransits 
(Lopez-Morales et al. 2006; Shan kland et al.ll2006[ ) which 
have thus far been unsuccessful. 

This paper discusses the effect of orbital parameters on 
the geometric transit probability of planets. We calculate 
orbital constraints that may be applied, particularly in 
the absence of transit signatures in photometric follow- 
up observations. Section 2 describes how the eccentricity 
and argument of periastron of known planetary orbits af- 
fect transit probability. It further presents applications 
of this effect to known RV planets and discusses how 
uncertainties in the orbital parameter values affect the 
reliability of the ephemeris calculations. In Section 3, 
we show how orbital constraints can be applied in the 
absence of a photometrically detected transit signal, and 
we discuss the potential transit yield and statistical sig- 
nificance of a scenario in which no transits are found in a 
large sample of RV planets. We summarize and conclude 
in Section 4. 

2. TRANSIT PROBABILITY 



Recent work by iBarned (|2007t ) and iBurkd 
showed that higher eccentricities of planetary orbits will 
increase their transit probabilities and, consequently, ex- 
pected yield for transit surveys. In this Section, we 
demonstrate the combined effect of the eccentricity and 
argument of periastron on transit probability. For ex- 
planations of the orbital parameters, includin g th e argu - 
me nt of per i astron w, we refer the reader to Kand (|2007l ) 
and iBarnei (j2007f ). We first explicitly derive the depen- 
dence of transit probability Pt as a function of eccen- 
tricity e, argument of periastron u, and orbital semi- 
major axis a (i.e., period). We discuss this dependence 
of Pt specifically with respect to u and period, apply 
the results t o a sa mple of 203 exoplanets compiled in 
iButler et ahl ([2006'), and briefly discuss how ephemeris 
calculations (and thus planning of photometric follow-up 
observations) are affected by uncertainties in e and uj. 
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2.1. Orbital Configuration 

For a circular orbit the geometric transit probabihty 
is proportional to the inverse of the semi-major axis, a, 
such that the inclination of the planet's orbital plane i 
must satisfy 



a cos i < Rp + 



(1) 



where Rp and Rj, are the radii of the planet and star re- 
spectively (jBorucki fc Summerslll984j ). For an eccentric 
orbit, the transit probability, Pj, may be expressed as 



Pt 



Rp + R* 
a(l — e cos E) 



(2) 



where e is the eccentricity of the orbit and E is the ec- 
centric anomaly. The eccentric anomaly and the true 
anomaly, /, are related to each other by 



cosE = 



cos/ 



1 + e cos / 



(3) 



where the true anomaly is defined as the angle between 
the direction of periapsis and the current position of the 
planet in the orbit. Equation [2] can then be evaluated 
at each point in the planetary orbit. The transit proba- 
bility can also be described in terms of the geometry of 
an ellipse. For an elliptical orbit, the separation of the 
planet and star is 

r^ (4) 
1 + ecos/ ^ ' 

As shown by iKand (|2007D . the place in a planetary or- 
bit where it is possible for a transit to occur (where 
the planet passes the star-observer plane that is per- 
pendicular to the plane of the planetary orbit) is when 
CO + f — tt/2. The transit probability may then be re- 
expressed as 



{Rp + -I- e cos(7r/2 - u)) 



(5) 



a(l-e2) 

consistent with the findings of lBarnei (|2007f ) . Equations 
[^landOboth yield the same result based upon the orbital 
configuration, but Equation [5] clearly shows the major 
role played by the the values of e and uj in determining 
the likelihood of a planet transiting the parent star. 

2.2. Argument of Periastron Dependence 

Equation[5]states the dependence of transit probability 
on the argument of periastron. As we rotate the semi- 
major axis of the orbit around the star we can observe 
how the transit probability varies. This dependence is 
shown in Figure [T] for eccentricities of 0.3 (dashed line) 
and 0.6 (dotted line) in comparison with the constant 
transit probability for a circular orbit (solid line). Since 
the shape of this variation is independent of period, P, 
the y-axes are scaled for both a 4.0 day and 50.0 day 
period orbits. Figure [1] assumes a Jupiter radius and 
a solar radius for the values of Rp and R^, respectively. 
Note that Pt scales linearly with the sum of these values 
(Equation El). 

The peak transit probability occurs at a; = 7i'/2, and 
the corresponding increase in Pt as compared to a circu- 
lar orbit can be significant: a factor of 1.5 for e = 0.3 
and a factor of 2.5 for e = 0.6. Moreover, the fraction of 
the orbital path which produces a higher value of Pt than 
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Fig. 1. — Dependence of geometric transit probability on the 
argument of periastron, uj, for eccentricities of 0.0, 0.3, and 0.6, 
plotted for periods of 4.0 days (left ordinate) and 50.0 days (right 
ordinate). Stellar and planetary radii are assu med t o be a Jupiter 
and solar radius, respectively. For details, see i|2.1l and. in partic- 
ular. Equation [5] 

the circular orbit with the same period (corresponding to 
the fraction of range in uj for which the dotted or dashed 
line is above the solid line in Figure [T]) increases with 
increasing eccentricity. 

The fraction of orbital orientations with e ^ pro- 
ducing lower transit probabilities than the correspond- 
ing circular orbits is made clear in Figure [2] in which a 
view from above the orbit pole of two planetary orbits 
is depicted. The range of uj in Figure [1] that produces 
lower values of Pt than a circular orbit corresponds to 
the angle between the intersection points shown in Fig- 
ure [2] for which the planet is located outside the circular 
orbit. For an eccentricity of 0.6 this angle is 9 = 105° 
and decreases with increasing eccentricity. However, the 
Keplerian nature of the orbit is such that, although the 
larger fraction of the orbital path is spent close to the 
star, the larg er fraction of time is spent farther away 
from the star (jBarne i|2003). This is a crucial aspect in 
designing a photometric follow-up campaign to monitor 
RV planets in eccentric orbits for possible transits. 

2.3. Period Dependence 

As demonstrated in Figure [TJ the peak transit prob- 
ability increases with eccentricity. Although the shape 
of Pt = /(t^) is independent of period, the magnitude 
of Pt changes as a function of period (Figure [J). Con- 
sequently, the fractional increase in Pt for eccentric or- 
bits can be su bstantial, as shown in ij2.2l and argued by 
lBarne£l(|2007l) . 

The current distribution of eccentricities for the known 
extra-solar planets indicates that orbits within 0.1 AU 
tend to be forced into nearly circular orbits through tidal 
circularization, whereas longer period orbits can po ssess 
a great range of eccentricities (jFord fc: Rasiol l2008h . In- 
deed most of the planets beyond 0.1 AU have eccentrici- 
ties in excess of 0.3. Thus, it is the longer-period planets 
whose transit probabilities are more likely to be affected 
by eccentricities than the short-period ones. 

In Figure [3] we show mean transit probability as a func- 
tion of period after averaging over < w < 27r, for the 
period range 1 < P < 50 days. Eccentricities of 0.0, 0.3, 
and 0.6 are shown with solid, dashed, and dotted lines, 
respectively. As expected, we see that doubling the ec- 
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Fig. 2. — A view from above the orbit pole of a circular (solid 
line) and eccentric (e = 0.6; dotted line) planetary orbit for = 
3tt/2. The angle 6 corresponds to the range of orbital orientations 
for which an elliptical orbit has a lower transit probability than a 
circular orbit with the same period (see ^2. H and Figure [Til. 
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Fig. 3. — The mean transit probability on a logarithmic scale, 
averaged over all values of lu (cf. Figure[T]l, as a function of period, 
for eccentricities of 0.0, 0.3, and 0.6. 

centricity from 0.3 to 0.6 creates a significant increase in 
the mean transit probability. Most affected are the longer 
period planets whose eccentric orbits can raise their like- 
lihood of transit from a negligible value to a statistically 
viable number for photometric follow-up. 

2.4. Application to Known Exoplanets 

If we assume circular orbits for each of the known ex- 
oplanets, the transit probability at intermediate to long 
period orbits makes photometric searches for planets in 
those regimes impractical. However, applying the orbital 
parameters of e and uj should in general lead to an overall 
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Fig. 4. — The geometric transit probability for a circular orbit 
(soli d curve ) along with the transit probability for 203 RV planets 
from ' Butler et a l. (2006) calculated from their orbital parameters 
(open circles). HD 17156b is indicated by a 5-pointed star. The 
sub-panel plots the difference in Pt between the the actual orbit 
and a hypothetical circular orbit for each of the planets. 

more favourable situation for transit detection. Depend- 
ing on the brightness of the host star and the cadence of 
the RV observations, a reasonable estimate of these two 
parameters is normally extracted from the RV fitting. 

Figure [3] shows the transit pro bability calcu l ated f rom 
orbital parameters provided by iButler et al.l ()2006[ ) for 
planets with estimates of e and u> (203 planets in total) . 
The transit probabilities are plotted against period, but 
are calculated from the semi-major axis, a, using Equa- 
tion O For the purposes of providing an approximate 
comparison of the relative transit probabilities, we as- 
sume a Jupiter and Solar radius for the values of Rp and 
Ri,, respectively. Hence, we can include the transit prob- 
ability for a circular orbit, shown in the figure as a solid 
line. In addition, the sub-panel in the plot shows the dif- 
ference in Pt between the actual orbit and a hypothetical 
circular one of the same period (residuals). The mean 
value of the residuals for all 203 planets is positive but 
relatively small (4.13 x 10~^), and is dominated by the 
low transit probability of the long period planets. The 
mean residual of planets with P < 100.0 days, however, 
yields an overall increase of ^ 0.5% in Pt. 

HD 17156b, a transiting planet with 21.2 day period 
(jBarbieri et al.l 2007). is shown as a 5-pointed star. Its 
transit probability is greatly increased by its orbital pa- 
rameters. Note that the actual Pt of HD 17156b is larger 
than the 5% shown in Figure [5] since the radius of the 
host star is 1.47 Rq. At longer periods, the planets with 
the largest residuals are HD 156846b, HD 4113b, and 
HD 20782b, which have periods of 359.51, 526.62, and 
585.86 days, respectively. The probability residuals for 
these three planets are 0.024, 0.032, and 0.025 respec- 
tively, the effect of which is to raise their transit prob- 
abilities to the same level as HD 17156b if it were in a 
circular orbit. It is worth noting that these three planets 
all have eccentricities close to 0.9 which is undoubtedly 
the primary cause of the increased transit probability. 

The increased transit probabilities of eccentric planets 
motivate photometric follow-up programs of RV plan- 
ets. Compared to transit surveys, such programs require 
much less telescope time since the time of transi t is, in 
principle, known. However, it was shown bv lKand (|2007f ) 
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that reliable constraints on e and w are needed to avoid 
significant offsets in predicted transit times. This is par- 
ticularly true of long period planets. In the case of plan- 
ets HD 156846b, HP 4113b, a nd HD 20782b, the un- 
certainties cited in But ler et al.l (2006) indicate that the 
values of e are all constrained to ±0.03 and the values of 
uj are constrained to ±3.0° (compared to the mean and 
median values for all 203 planets of 6oj = 20° and 10°, 
respectively). Thus, ephemerides for these planets can 
relatively reliably be determined from RV fit parameters 
alone. The transit duration is on the order of 12 hours 
for these planets, ensuring that one will practically never 
observe both the ingress and the egress of the transit dur- 
ing a single orbit from the ground. However, the large 
transit duration and relative low uncertainties in e and 
Lo will increase the chances of observing at least a partial 
transit during the predicted observing window. 

3. CONSTRAINING ORBITAL PARAMETERS 

In §2, we discussed transit probability as a function of 
various system parameters as well as aspects of poten- 
tial photometry follow-up campaigns. Here we focus on 
what can be learned from the presence and absence of a 
planetary transit in follow-up observations. 

For a transiting planet the physical properties (such as 
the mass, radius, and density) can be calculated, leading 
to determination of (as opposed to constraints on) system 
parameters of the planet. Furthermore, the orbital incli- 
nation can be corn pared with the plane of stellar rotation 
(|Winn et al.ll2007D and used to test planetary models re- 
garding co-planar orbits. 

However, even the absence of a planetary transit sig- 
nature in photometric data can lead to interesting con- 
straints on the orbital parameters. Below, we elaborate 
on these c onstraints and appl y the results to the afore- 
mentioned |Butler]e£^ ()2006[ ) sample of RV planets. 

3.1. Orbital Radius versus Stellar Radius 

One implicit as sumption in the derivation of the tran- 
sit probability bv lBorucki fc Summers! () 19841 ) is that the 
planet remains well outside the star in order to produce 
the s olid angle of the planet's shadow (see also iBarnesI 
l2007f) . As a result, the calculation of Pt in Equation[5]be- 
comes invalid for extreme orbits with a small semi-major 
axis and a high value of eccentricity (see H2.1^ . 

To quantify this assumption, we use Equation |4] to cal- 
culate the maximum eccentricity, emax, allowed as a func- 
tion of the planet-star separation in units of a/R where 
R = _Rp + i?^. Applying the constraint r > R when 
f = (i.e., the planet is outside the star at periapsis) to 
Equation [5] results in 



a(l 



and thus. 



1-^. 

a 



(6) 
(7) 



Equation [7] is plotted in Figure [5] for values of a/R 
ranging from 1 to 30. Also shown are dot-dashed 
lines which indicate the a/R v alues fo r OGLE-TR-56b 
(jKonacki et all [Ml . X0-5b (jBurkd r2008), and HD 
17156b. 

The restrictions upon the maximum eccentricity begin 
to become significant for a/R < 10, which encompasses 
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Fig. 5. — The maximum orbital eccentricity, emax, plotted as a 
function of the planet-star separation in units of a/R (see Equa- 
tion [7]l in order for a planet to remain outside the surface of its 
parent star. This requirement is purely based on system geometry 
and does not take into account tidal effects or planet-planet inter- 
actions, but requires that a > R when / = 0. Dot-dashed lines 
indicate values for a/R for various transiting planets. 

most of the known transiting exoplanets. This restriction 
is purely based upon orbital dynamics and there are un- 
doubtedly additional limitations on the eccentricity due 
to tidal effects in this region. 

3.2. Orbital Inclination and Argument of Periastron 

One of the primary advantages of observing an exo- 
planet transiting the host star is that it eliminates the 
ambiguity in the planetary mass created by the unknown 
orbital inclination angle, i. The precise value of the in- 
clination can be derived from the impact parameter of 
the transit across the stellar disk, defined by 



b = 



a cos I 
R. 



(8) 



and measurable from th e shape of the lightcurve and th e 
planet-star radius ratio (Sea ger fc Mallen-Ornelaj|2003f) . 
Due to the constraint placed on i by the presence of tran- 
sits, the true planetary mass will be within a few percent 
of the value as determined from RV measurements alone. 

The data available for transiting planets from 
the Extra-solar Planet s Encyclopaedia^ and from 
[Torres. Winn, fc Holmai] ([2008) show that the current 
distribution of inclination angles extends from 90°to al- 
most 78°. The transiting planets whose orbits feature 
numerically lower values of i (i.e., more "face-on") are 
dominated by the very hot Jupiters, such as OGLE-TR- 
56b w hich has an inclination of 78.8° ± 0.5° (jPont et al.l 
l2007l) . 

If, however, a planet is determined not to transit, then 
limits may be placed upon the orbital inclination if the 
eccentricity and argument of periastron are known from 
RV measurements. This results from re-expressing Equa- 
tion [2] as follows: 



cos 2 



Rp + R* 
a{l — ecosE) 



(9) 



Figure [S] shows the maximum inclination for various val- 
ues of period, e, and to. These are calculated by holding 

http://cxoplanct.eu/ 



Constraining Orbital Parameters 



5 



e = 0,00 
e = 0.30 
e = 0,60 



lO II 



0.5 



1 



1.5 



Argument of periastron (tt radians) 

Fig. 6. — The maximum orbital inclination for a non-transiting 
planet as a function of the argument of periastron, w, for eccen- 
tricities of 0.0 (solid line), 0.3 (dashed line), and 0.6 (dotted line), 
plotted for periods of both 4.0 days and 50.0 days. 

period and e fixed whilst varying uj using Equations [3] 
and [9l We further assume a Jupiter radius and a solar 
radius for the values of Rp and i?* respectively. For non- 
transiting planets on orbits with e ^ whose periastron 
is aligned towards the observer (i.e., lj ^ 7r/2), Figure 
[6] shows that the constraint on the inclination can be as 
high as i < 75°, depending upon the orbital period. This 
is particularly useful for those planets whose mass esti- 
mate places them close to the brown dwarf regime. Note 
that Equation [H] reduces to Equation [7] when i — 90° and 
LO = tt/2, consistent with the requirement that the planet 
remain outside the star. 

3.3. Orbital Inclination and Eccentricity 

As we show in §3.21 the fact that a planet is found to 
not transit limits the possible combinations of e, w, and 
i. We now consider what constraints may be placed upon 
the orbital inclination for a non-transiting RV planet as a 
function of eccentricity for the specific examples of when 
the periapsis is aligned towards (w ~ 7r/2) and away from 
(w ~ 37r/2) the observer. 

Figure [7] illustrates the range of orbital inclinations 
that are excluded for two orbits (shown edge-on) of 
non-transiting planets. Both of these orbits have the 
same semi-major axes but different eccentricities and are 
aligned such that w = 37r/2 (i.e., periapsis occurs behind 
the star as seen from the observer). In this case, the 
range of possible values of i increases with decreasing or- 
bital eccentricity ((/>! > ^2)- The opposite is true when 
w = 7r/2. In fact, the inclination in that case is only 
constrained by the requirement that the planet remain 
outside the star during periapsis (Equation [7]). 

Figure [8] graphically demonstrates these constraints by 
plotting Equation [51 except now we fix the period and 
UJ and vary e. A Jupiter radius and a solar radius are 
assumed for the values of Rp and i?* respectively. The 
lines in these plots represent the maximum values for 
i for a non-transiting planet as a function of e. These 
calculations are performed for four different periods and 
the two aforementioned orientations of w: u; — tt/2; left 
panel; case (a), and to = 37r/2; right panel; case (b). It 
is worth noting that case (a) in Figure [8] represents the 
physical constraint described in Figure [5l since Equation 
[HI reduces to Equation [7] for uj — tt/2, as stated above. 

In case (a), for example, a non-transiting planet in a 
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Fig. 7. — An edge-on view of two planetary orbits with the 
same values of semi-major axis, showing the range of excluded 
inclinations for an orbit with low eccentricity (dotted lines; tpi) 
and high eccentricity (dashed lines; 02 ) for which the planet does 
not transit the parent star. For lj = 37r/2, the range of possible 
inclinations decreases with increasing eccentricity. 

4-day orbit with e = 0.4 has a range of possible incli- 

o 

nation angles of z < 80 . For larger values of e (at any 
period), the periastron distance of the planet will be- 
come so small that almost all values of i are possible, 
the maximum value of e at each period being defined by 
Equation [T] The dependence of i upon e is weaker for 
case (b) (note different scale for the left panel in Figure 
[5]) since the periastron passage now happens behind the 
star (Figured]) as seen from the observer, and thus, the 
range of possible i-values is not very constrained by e. 

We now consider a planet discovered using the transit 
method with known i, a, Rp, and i?*, but unknown values 
for e and uj. Is it possible to constrain e in this scenario? 
Using Equation [9l the eccentricity can be expressed as 
follows 



1 



cosE 



1 



Rp + R* 



a cos I 



(10) 



However, there exists a degeneracy between e and w = 
f{E) such that one cannot place constraints on one pa- 
rameter without knowledge of the other. Additionally, as 
shown in i )3.21 the constraint upon i is only limited by the 
orbital boundary defined by Equation [7] when uj ~ tt/2. 

Therefore, a meaningful constraint may only be placed 
upon e for values of uj for which the orbital inclination is 
greater than the maximum predicted for a circular orbit 
(i.e., the region above the solid line shown in Figure O. 
For case (b) [uj — 37r/2), Equation [TPl reduces to 



Rn + R* 



- 1. 



(11) 



As an example, consider the two known transiting plan- 
ets TrES-3 and TrES-4. The fit parameters shown in Ta- 
ble [T] for the values of i?^, Rp, i, and a are those reported 
by th e discovery papers for TrES-3 (lO 'Donovan et aLl 
I2007f) and TrES-4 (jMandushev et al.ll2007D . Also shown 
in Table [1] are the maximum eccentricities for both case 
(a) and case (b) . For case (a) , the maximum eccentricity 
is ~ 0.8 for both planets. For case (b), the maximum ec- 
centricity for these two planets is 0.4-0.5 and are plotted 
in the right panel of Figure \8\ In each case, the max- 
imum eccentricities arc remarkably similar because the 
longer period of TrES-4 is compensated by the relatively 
large radii of the star and planet. 

3.4. Global Statistics 

The total number of transiting planets discovered thus 
far via radial velocity surveys docs not necessarily reflect 
the true number of transiting planets in this sample. At 



6 Stephen R. Kane & Kaspar von Braun 




Fig. 8. — Maximum orbital inclination as a function of e for non-transiting planets, plotted for four differen t per iods (see Equation |9]l . 
The left panel, case (a), is for oj = 7r/2 and is based on the requirement that the planet remain outside the star ( i|3 11 1. The right panel, case 
(b), shows the situation for a; = 37r/2, and is based on the geometrical arguments outlined in i|3.3l and shown in Figure [7] The location of 
the maximum eccentricities of the known transiting planets TrES-3 and TrES-4, given their parameters in Table[l]and assuming uj = 3it/2, 
are indicated by 5-pointed stars. 

TABLE 1 ° 

Fit parameters for TrE S-3 and TrES-4 fro m 

lO 'Donovan et al] (I2007D and IMandushev et al] Il2007r 
respectively, along with the calculated maximum s . 

eccentricities. ° ■ 



Parameter 



TrES-3 



TrES-4 



Hp 
i 

a 

e{uj = 7r/2) 
e(i.j = 37r/2) 0.404 ± 0.041 



0.802 ± 0.046 
1.295 ±0.081 
82.15 ±0.21 
0.0226 ± 0.0013 
0.808 ± 0.083 



1.738 ±0.092 
1.674 ±0.094 
82.81 ±0.33 
0.0488 ± 0.0022 
0.818 ±0.073 
0.451 ±0.040 



the time of writing, most of the known radial velocity 
planets have not been adequately monitored photomet- 
rically in order to rule out transits. We can estimate the 
number of planets that should be transiting, and deter- 
mine the significance of a hypothetical null result from 
a photometric follow-up campaign by applying the re- 
sults of this paper to the iButler et al.l (p006 ) RV planet 
sample. 

The host star prope rties and planetary orbital param- 
eters provided by Butler et al.l (j2006) form the founda- 
tion of a Monte-Carlo simulation of the transit proba- 
bilities calculated from Equation [5l The planetary radii 
Rp are assumed to be one Jupiter radius as used in pre- 
vious sections. However, the stellar radii i?^ are esti- 
mated indi vidually from the values oi B — V provided by 
IButler et a l. (2006), assuming the host stars are dwarf 
stars ()Co^ ,2000,1 . The or bital elemen t s a, e and w 
are directly extracted from IButler et al.l (|2006D . Using 
these values, we calculate Pt for each of the 203 stars in 
the sample and randomly determine if the planet tran- 
sits. This yields an integer number of projected tran- 
sits from the sample. By performing these calculations 
~ 100000 times, we produce a probability distribution 
for the number of transiting planets expected from this 
sample, shown in Figure (9] 

The simulated probability distribution has a mean 
value of ~ 4.5 transits peaking at Pt ~ 0.2 with a stan- 
dard deviation of ^ 2.0. For comparison, we also gener- 




Number of transiting planets 

Fig. 9. — T he probability distri bution (solid hne) for the 203 
planets in the IButler et al.l I I2006I) sample, predicting the num- 
ber of transiting planets based on their estimated orbital param- 
eters. Over-plotted is a gaussian distribution (dashed line) using 
the mean and standard deviation of the simulation results. 

ated a gaussian distribution profile using this mean and 
standard deviation. The a priori probability that none 
of the planets in this sample transit their host star is 
~ 1%. In fact, three of the planets in this sample are 
known to transit, specifically HD 17156b, GJ 436b, and 
HD 147506b. Hence the current number of transiting 
planets from this sample is almost Icr below the expec- 
tation. 

We further note that the sample of RV planets is biased 
toward numerically higher values of i since detection effi- 
ciency will increase with higher i for a given RV precision. 
As such, the expected number of transiting planets in 
the sample should be regarded as a lower limit. Though 
the discrepancy between known and expected transiting 
planets is not significant in this low-number regime, it is 
nevertheless quantifiable, and we conclude that further 
transit discoveries in this sample are possible or even 
likely. Any such additional detections would, in turn, 
lead to further understanding of the respective observa- 
tional biases of the RV and transit methods. For exam- 
ple, the observational bias leads to an observed difference 
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between the period distributions of planets discovered by 
the transit method a nd the radial velocity method, as 
discussed in detail by iGaudi. Seager. &: Mallen-OrnelasI 
(|2005i) . 

4. CONCLUSIONS 

It is still uncertain at this stage how many of the known 
radial velocity planets transit their parent stars. What 
is clear is that the eccentricity distribution of the known 
exoplanets will increase the transit likelihood, making 
detections for long-period planets, such as HD 17156b, 
feasible. We have shown in this paper that there is 
enough potential amongst longer period planets for tran- 
sit detections to motivate a photometric monitoring cam- 
paign at the p redicted times of transit for these targets. 
iFleming et al.l (j2008.) have shown that long-period tran- 
siting planets may yet be discovered through ground- 
based transit surveys, particularly if data sets from dif- 
ferent surveys are c ombined . 

As pointed out by Barnes' (2007), eccentric planets that 
have a periastron oriented away from the observer are far 
more likely to exhibit a secondary than a primary eclipse. 
The detection of such a secondary eclipse is considerably 
more challenging than for a primary eclipse since it re- 
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lies on a minimum level of planetary flux and is best 
pursued at infrared wavelengths. The discussion in §3.31 
shows that even an assumption of a; = 37r/2 can place 
constraints on the orbital inclinat ion. A prime can didate 
for such a study is HD 80606b (Nae f et al.ll2001h which 
has a period of 111.87 days and an eccentricity of 0.927. 
Scaling Figured] to this period and eccentricity yields a 
secondary transit probability of ~ 15%. 

Many of the results presented in this paper can eas- 
ily be applied to any system since the results generally 
scale linearly with the sum of the stellar and planetary 
radii. Through applying these results to current and fu- 
ture radial velocity planet discoveries, one can choose 
targets for an efhcient observing campaign which may 
help to discover long-period transiting planets and hence 
add invaluable information to planetary structure and 
formation theories. 
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